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Abstract 

In this paper, we introduce a new class of distributions by compounding the exponen¬ 
tiated extended Weibull family and power series family. This distribution contains sev¬ 
eral lifetime models such as the complementary extended Weibull-power series, generalized 
exponential-power series, generalized linear failure rate-power series, exponentiated Weibull- 
power series, generalized modified Weibull-power series, generalized Gompertz-power series 
and exponentiated extended Weibull distributions as special cases. We obtain several prop¬ 
erties of this new class of distributions such as Shannon entropy, mean residual life, hazard 
rate function, quantiles and moments. The maximum likelihood estimation procedure via 
a EM-algorithm is presented. 

Keywords: EM-algorithm, Exponentiated family, Maximum likelihood estimation, Power se¬ 
ries distributions. 


1 Introduction 


The extended Weibull (EW) family contains variou s well-kiiown distributions such a s exponen¬ 


tial, Pareto, Gompertz, Weibull, linea r failu r e rate (IBarlowl. 


2003|), additive Weibull ( Xie and Lai ^ 


distributions. For more details see 
Using the given method by 


1995 


Almalki and Yuan. 


Nadaraiah and Kot! 


3 ( knn.jf 


1968), mo dified Weibu’ 


201^1 and Chen (IC 


I (ILai et al. 


I) and IPha.m and I^ail ( 20071) . 


nen 


2OOOI) 


Gupta and Kundul (119991) . the EW family can be generalized. 


We call it exponentiated extended Weibull (EEW) distribution. The cumulative distribution 
function (cdf) of this distribution is 


G(x; a,/3, 0 ) = [1 — e a > 0, /3 > 0, x > 0, 


( 1 . 1 ) 
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and its probability density function (pdf) is 


g{x- a, /3, 0) = a/3/i(x; [1 - 


( 1 . 2 ) 


where 0 is a vector of parameters, and H(x-,Q) is a non-negative, continuous, monotone 
increasing, differentiable function of x such that H{x]&) —0 as x ^ O"*" and H{x-, 0) 
as X ^ oo. It is denoted by EEW(q!, /3, 0). 


oo 


The EEW distribuyion is a flexible family anc 


modified Weibull (Carrasco et al 


conta ins many exponentiated distributions 


such as generalized exponenti al (IGuota and Kundu 
I 993 I I. generalized Rayleigh (ISurles and Padeettl . 


2001 


1999I'). exponentiated Weibull (jMudholkar and Srivastaval . 


Kundu and Ragatj 


2008I1. generalized lin ear failure rate (jSarhan and Kundu . 


20091 ). and generalized Gompertz (jEl-Coharv et al.l . l201.'ll ) distributions. 


2 OO 5 II. generalizec 


In recent years, many distributions to model lifetime data have been introduced. The ba¬ 
sic idea of introducing these models is that a lifetime of a system with N (discrete random 
variable) components and the positive continuous random variable, say Xi (the lifetime of ith 
omponent), can be denoted by the non-negative random variable Y = min(Xi,..., Xjv) or 
Y = max(Xi,..., X^r), based on whether the components are series or parallel. 

In this paper, we compound the EEW family and power series distributions, and introduce 
a new class of distribution. This class of distributions can be applied to reliability prob¬ 
lems and its some properties are investigated in this paper. We call it exponentiated extended 
Weibull-power series (EEWPS) class of distributions. In similar way, so me distributions are pro¬ 
posed in literature: The exponential-power series (EP) distribution by 


(|200^ b Weibull-power series (WPS) distributions by 
eralized exponential-power series (GEP) distribution W 


Morais and Barreto-Souz. 


plementary exponential power serie s 


(EWPS) distribution by 


Silva et al 


Bidram and Nekoukhoul (|2013l l. Burr-power series by 




kL 


Elores et al 


(|2ni^ ). double b ounded Kumaraswamy-po wer series by 


Mahmoudi and Jafari 


Ghahkandi and Ganiali 

T 


201 ih . gen- 
(]2012 l. com- 


([20131), extended Weibull-power series 


linear failure rate-power series (GLER P) distribution by 


Saunde rs-power series distribution b y 


ries by 


Bourguignon et al 


Silva and Cordeirol (l2()13lh g eneralized 
Alamatsaz and ShamJ ( 2ni4 ). Birnbaum- 


(j2ni4l ). linear failure rate-power se- 


Mahmoudi and Jafaril (j2014l l. and complementary extended Wei 


Gordeiro and Silva 

(2014 

). Similar procedures are 

(2011 

), 

Nadaraiah et al. 

(2014a 

) and 

Louzada et al. 


3 ull-po w er series by 


Roman et al 


tributions with discrete distributions, 


11201 


Nadaraiah et al 


pounding in R software (jP, Development Gore Team 


(l2mi). 


Lu and Shi 


20141 1 


I). Eo r compounding continuous dis- 


|2013l l introduced the package Com- 































































































































We provide three motivations for the EEWPS class of distributions, which can be ap¬ 
plied in some interesting situations as follows: (i) This new class of distributions due to the 
stochastic representation V = max(Xi, ... ,Xjy), can arises in parallel systems with identical 
components, where each component has the EEW distribution lifetime. This model appears in 
many industrial applications and biological organisms which the lifetime of the event is only 
the maximum ordered lifetime value among all causes, (ii) The EEWPS class of distributions 
gives a reasonable parametric fit to some modeling phenomenon with non-monotone hazard 
rates such as the bathtub-shaped, unimodal and increasing-decreasing-increasing hazard rates, 
which are common in reliability and biological studies, (hi) The time to the last failure can be 
appropriately modeled by the EEWPS class of distributions. 

The remainder of this paper is organized as follows; The pdf and failure rate function of the 
new class of distributions are given in Section [21 The special cases of the EEWPS distribution 
are considered in Section [3l Some properties such as quantiles, moments, order statistics, 
Shannon entropy and mean residual life are given in Section [H Estimation of parameters by 
maximum likelihood are discussed in Section |5l Application to a real data set is presented in 
Section El 


2 Introducing new family 


A discrete random variable, is a member of power series distributions (truncated at zero) if 
its probability mass function (pmf) is given by 


Pn = PiN = n) = 


anX^ 

C(A) 


, n = 1,2 ,..., 


( 2 . 1 ) 


where an > 0, C'(A) = ^ anX^, and A G (0, s) is chosen in a way such that C'(A) is finite 


n=l 


and its first, second and third derivatives are defined and shown by C'(. ). C"(. ) and C'"(.) 


respectively. The term “power series distribution” is generally credited to 


NoackI ( 195(1 ). This 


family of distributions includes many of the most common distributions, including the binomial, 
Poisson, geometric, neg ative binomial, logari thmic distributions. For more details about power 


series distributions, see 


Johnson et al 


((20051), page 75. 


Theorem 2.1. Let N be a random variable denoting the number of failure causes which it is a 
member of power series distributions with pmf in m- Also, For given N, let Xi,X 2 , ■.■,Xn be 
independent identically distributed random variables from EEW distribution with pdf in dLH). 
Then = maxi^j^ArlW} has EEWPS class of distributions is denoted by EEWPS (a,/3, A, 










0) and has the following pdf: 


fix) = apXhix-, 

Proof. The conditional cdf of -^(at) | = n has EEW(a, n/3, 0). Hence, 

n 

P(X(^) <x,N = n) = -^[l- 
and the marginal cdf of ^(at) is 

, C(A(1 

=-C(A)-• " > “■ 


X > 0. (2.2) 


(2.3) 


(2.4) 


The derivative of F with respect to x is (|2.2p . Therefore, has EEWPS distribution. □ 

Proposition 1. The pdf of EEWPS class can be expressed as infinite linear combination of 
density of order distribution, i.e. it can be written as 


fix) = 'P^Png{n)ix;a,nj3,&), 


n=l 


where 5 („)(x; a, nj3, 0) is the pdf of EEW distribution with parameters a, n/3 and 0. 
Proof. Consider t = 1 — go 

fix) = a(3\hix- ciX) 

OO 

Y, nan(At^)"“^ 


C(A) 


anA" 


E ^na/3h(x; ©)e-“^^(-;®)r/3-i 


n 

OO 


'P^Png{n)ix\wnld,&). 


n=l 


Proposition 2. The limiting distribution o/EEWPS(/3, A, 0) when A —)• 0"*“ is 


(2.5) 


□ 


lim Fix) = [1 - 
A^0+ 


which is a EEW distribution with parameters a, c(3 and 0, where c = min{n G N : On > 0}. 








Proof. Consider t = 1 — e So 


lim F(x) 
A^o+ 


1 - C{Xt0) 
hm 

A— >-0+ C(A) 


lim 

A^0+ 


CXD 

E 

n=l 


oo 


E anX^ 

n=l 


= lim 
A->'0+ 


OO 

acF^+ E 

n=c-\-l 

oo 

flc + E anX^-^ 

n=c+l 


= P^. 


□ 


Proposition 3. The hazard rate function of the EEWPS class of distributions is given by 


aX(3h{x]&){l — t)t^ ^C'[Xt^) 
C(A) - C{Xt^) 


( 2 . 6 ) 


where t = 1 — e 


Proof. Using (|2.2I) . (|2.4I) and definition of hazard rate function as r{x) = /(x)/(l — F{x), the 
proof is obvious. □ 


3 Special cases 


In this Section, we consider some special cases of the EEWPS distribution. 


3.1 Complementary extended Weibull power series 

If /? = 1, then the pdf in (|2.2I) becomes to 




(3.1) 


which is the pdf of complementary extended Weibull power series (CEWPS) class of distribu¬ 


tions introduced by 


Cordeiro and Silval (120141 1. 


3.2 Generalized exponential-power series 

If F[{x] 0) = X, then the pdf in (|2.2I) becomes to 


C' (X(l — 

/(x) = a/3Ae-"^(I - ^ ^ x > 0. 

t;(A) 


(3.2) 


which is the pdf of generalized expon ential-power series (GEPS) class of distributions intro¬ 


duced by 


Mahmoudi and Jafaril (j2012l b The GEPS class of distributio ns contains com p 


tary exponentiated exponential-geometric distribution introduced by 


Louzada et al 


emen- 


(1201311 

















exponential distribution introduced by 


Cancho et al 


(120111 1 and 


complementary exponential-geometric d istribution introduce d bv iLouzada et alJ (1201111. Poisso n- 

i2Qllll 


Louzada-Neto et al 


compl ementary exponential -power series class of distri butions introduced 


(|2013l l. generalized exponential distribution introduced by^ 
eralized exponential-geometric distribution introduced by 


Gupta and Kundu 


Bidram et al. 


(1201011 . 




Flores et al 


I 999 I I and gen- 


3.3 Generalized linear failure rate-power series 

If = ^ + then the pdf in (|2.2p becomes to 

C (A(l-e- 


f{x) = j3X{a + hx)e 


■(1-e 


—ax—1 


'-f 


C{\) 


X > 0. 


(3.3) 


which is the p df of generalized linear f ailure rate-power series (GLFRPS) class of distributions 


introduced by 


Alamatsaz and Shami 


(I 2 OI 4 II. It is a mod ification of generalized linear failure 


rate distribution introduced bv ISarhan and Kundul (1200911 an d generalized linear failure rate- 


geometric distribution introduced by 


Nadarajah et al 


( 2014bl ). If 6 = 0, it becomes to GEPS 


cla ss of distributions. Also, If 0 = 1, it becomes to linear failure rate-power series introduced 


by 


Mahmoudi and Jafaril (j2014l l. 


3.4 Exponentiated Weibull-power series 

If H{x-, 0) = x'*', then the pdf in ()2.2I) becomes to 


fix) = a/3A7x^-'e-“""(l - ^ ^ > q. 


(3.4) 


wh ich is the pdf of exponen t iated Weibull-power series (EWPS) class of distributions introduced 


by 


introduced by 


Mahmoudi and Shiran! (1201^1. It is a modi fication of exponentiated Weibull distribution 


Mudholkar and Srivastaval (Il993ll. 


ometric distribution introduced by_ 
Weibull distribution introduced by 


Toieiro et al. 


t is co ntain the complementary Weibull ge- 
J 2014 b Also, the Marshall-Olkin extended 


Gordeiro and Lemontd (|2013ll is a special case of EWPS. 


3.5 Generalized modified Weibull-power series 

If FI(x; 0) = x'>'exp(Tx), then the pdf in (|2.2p becomes to 

fix) = a(3\x'^~^i'y -|- rx)e 


„P(„) C-(A(1 

(1 — e-"^'^®^P(™))i-(^C(A) ’ 


X > 0, (3.5) 







































































and we call generalized modified Weibull-power series (GMWPS) clas s of distributions. It is 


contained the generalized modified Weibull distribution introduced by 


Carrasco et al. 


( 1200811 . 


If T = 0, then GMWPS class of distributions becomes to EWPS class of distributions. 

3.6 Generalized Gompertz-power series 

If H{x] 0) = — I), then the pdf in (12.21) becomes to 

f{x) = aj3\e^^e ^^(1 —e ^x > 0. (3.6) 

C(A) 

and we call generalized Gompertz-power seri es class of distributions . It is contained the gen¬ 


eralized Gompertz distribution introduced by 


El-Goharv et al 


(I2OI3I I 


4 Statistical properties 

In this section, some properties of EEWPS class of distributions such as quantiles, moments, 
order statistics, Shannon entropy and mean residual life are derived. Using (|2.5I) . we can obtain 

00 00 

F{x) = ^pnG(„)(x;a,n/3,0) = , (4.1) 

n=l n=l 

where t = 1 — . Based on the mathematical quantities of the baseline pdf (x; a, n/3, 

0 ), we can obtain some statistical quantities such as ordinary and incomplete moments, gen¬ 
erating function and mean deviations of this family of distributions. 


4.1 Quantiles and Moments 


Let 


X = G 


-1 


C-MG(A)U) ^ 


where U has a uniform distribution on (0,1), G~^{y) = ln(l — y~^)] and C“^(.) is the 

inverse function of C(.). Then X has the EEWPS(q;, /?, A, 0) distribution. This result helps in 
simulating data from the EEWPS distribution with generating uniform distribution data. 


Theorem 4.1. Consider X ~ EEWPS(a,/3, A, 0). Then the moment generating function of 
EEWPS is 



(-l)^My(t), 


(4.2) 


where Y has EEW(a(j -|- 1), 1, 0). 










Proof. The Laplace transform of the EEWPS class can be expressed as 


L{s) = E(e-*^) = ^ P(iV = n)Ln{s), 

n=l 


where Ln{s) is the Laplace transform of EEW distribution with parameters a, n/? and 0 given 
as 


Lnis) 


/*+oo 

/ e-*^na/3/i(x; [1 - 

Jo 


na(3 

OO 

^n/3 








j + i 




where ^ 1 ( 5 ) is the Laplace transform of the EEW(q;(j + 1),1,0). Therefore, the moment 
generating function of EEWPS is 


Mx{t) 


00 

^ ^ PnL'n ( 0 

n=l 


n=l j=0 y 

n=l j=0 V / 


□ 


Theorem 4.2. The noncentral moment functions of EEWPS is 


“ 5 cw g 0+ ij =S g"”0+1 > 


(4.3) 


where y). = E\Y'^] and Y has EEW(a(j + 1), 1, 0). 


Proof. We can use Mx{t) to obtain /r^- But from the direct calculation, proof is obvious. □ 

With considering H{x) = x'^ and C'(A) = A(1 — A)“^, we calculated the first four moments 
with different values of parameters for the EEWPS distribution using ()4.3I) . Also, we computed 
these values from the direct definition by numerical integration. We found that the results are 
same. The values are given in Tables [H 




Table 1: The four moments of EEWPS mo del. 


a 

/? 

A 

7 

Ml 

M 2 

M 3 

M 4 

0.3 

0.3 

0.2 

2.0 

0.936 

1.594 

3.520 

9.164 

0.3 

0.3 

0.2 

5.0 

0.856 

0.884 

1.011 

1.237 

0.3 

0.3 

0.8 

2.0 

1.656 

3.719 

9.722 

28.292 

0.3 

0.3 

0.8 

5.0 

1.150 

1.446 

1.916 

2.631 

0.3 

2.0 

0.2 

2.0 

2.192 

5.446 

14.976 44.841 

0.3 

2.0 

0.2 

5.0 

1.345 

1.853 

2.606 

3.734 

0.3 

2.0 

0.8 

2.0 

2.835 

8.733 

28.694 99.531 

0.3 

2.0 

0.8 

5.0 

1.500 

2.285 

3.530 

5.521 

00 

o 

0.3 

0.2 

2.0 

0.573 

0.598 

0.808 

1.289 

00 

o 

0.3 

0.2 

5.0 

0.704 

0.597 

0.561 

0.565 

00 

o 

0.3 

0.8 

2.0 

1.014 

1.394 

2.232 

3.979 

00 

o 

0.3 

0.8 

5.0 

0.945 

0.977 

1.064 

1.201 

00 

o 

2.0 

0.2 

2.0 

1.342 

2.042 

3.439 

6.306 

00 

o 

2.0 

0.2 

5.0 

1.106 

1.252 

1.446 

1.704 

00 

o 

2.0 

0.8 

2.0 

1.736 

3.275 

6.589 

13.997 

00 

o 

2.0 

0.8 

5.0 

1.233 

1.543 

1.960 

2.519 

2.0 

0.3 

0.2 

2.0 

0.362 

0.239 

0.204 

0.206 

2.0 

0.3 

0.2 

5.0 

0.586 

0.414 

0.324 

0.271 

2.0 

0.3 

0.8 

2.0 

0.641 

0.558 

0.565 

0.637 

2.0 

0.3 

0.8 

5.0 

0.787 

0.677 

0.614 

0.577 

2.0 

2.0 

0.2 

2.0 

0.849 

0.817 

0.870 

1.009 

2.0 

2.0 

0.2 

5.0 

0.921 

0.867 

0.835 

0.819 

2.0 

2.0 

0.8 

2.0 

1.098 

1.310 

1.667 

2.239 

2.0 

2.0 

0.8 

5.0 

1.026 

1.070 

1.131 

1.210 


4.2 Order statistic 

Let Xi,X 2 ,..., Xm be a random sample of size m from EEWPS(a, /?, A, 0), then the pdf of 
the ith order statistic, say Xi^m, is given by 


fi-.mix) — 


m\ 




'c{xtf^y 

i-l p 

[ C{X) \ 



m\ 


E 


(i - l)!(m V J 


ml 


j=0 
oo m—i 


m — I 


1 - 


{-ly 




C{\) 


(i — l)!(m — i)! ^ ^ 

' ^ n=l 7=1 


X] E P^3{n) {x; a, nl3,&) 


ml 


j=0 
oo m—i 


m — I 


{-ly 


C{Xt^) 

L“^J 


j+i-l 


_ iVfm - Tv E E ^ jPn 9 in ) ix ; a , n / 3 , 0) 

V n=l j=0 


L“^ 


j+i-l 



where / is the pdf of EEWP class of distributions, t = 1 — e and wj = *)(—1)-^. 

Also, the cdf of Xi-m is given by 

m m—k / / \ 


k=i j=0 


m — k\ /m\ rC'(At^) 
j )[k) 


An analytical expression for rth moment of order statistics Xi-m is obtained as 

ml 


oo m—i 


E\XIJ = -T 


vEE w,PnE[Z^{F{Z)y+^-\ 


ii — l)!(m — i)! ^ ^ 

^ ^ '' n=l j=0 

where Z has a EEW distribution with parameters a, nj3 and 0. 


4.3 Shannon entropy and mean residual life 

The maximum entr opy met h od is a powerful technique in the field of probability and statistics. 


It is introduced by 


JavnesI (|l957l i and closely related to the Shannon’s entropy. Also, it is 


applied in a wide variety of fiel ds and used for the characterizatio n of pdf’s; see, for example , 
Kapui ( 1994 1 Soofi ( 2000 1 and Zografos and Balakrishnan ( 2009 1. Shore and Johnson ( 1980 1 
treated the maximum entropy method axiomatically. 

Considers a class of pdf’s 


F = {fix; a, /3, A, 0) : FfiT^X)) = I3i, i = 0,1,...., m}, 


(4.4) 


where Ti(A),..., Tm(Ar) are absolutely integrable functions with respect to /, and To(A") = 1. 
Also, consid er the shannon’s entropy of none-negative continuous random variable X with pdf 
/ defined by IShannonI (j 194811 as 


r+oo 

Hshif) = F[- log fix)] = - / fix) logif ix))dx. (4.5) 

Jo 

The maximum entropy distribution is the pdf of the class F, denoted by f^^ determined as 

/^®(x; A,/3,0) = aig max Hshif)- 

Now, suitable constraints are derived i n orde r to provide a maximum entropy character¬ 
ization for the class (|4.4|) based on Ijavn^ ( 1957 1. For this purpose, the next result plays an 
important role. 

Proposition 4. Let X has EEWPS(a,/3, A, 0) with the pdf given by (12.2p . Then, 

i. 

A 


E 


log(C''(A(l - e 


_ -aH{X-,& 




CiX) 


E 


C'(A(1 - 




































X log(C''(A(l — e 


-aH(Y-,@ 




ll. 


III. 


E[iog{h{x-,e))] = 


C{X) 


E 


E 




C(A) 


E 


C'(A(1 - log(/i(y; 0)) 


C'(A(1 - log(l - 


where Y follows the EEW distribution with the pdf in O)- 

An explicit expression of Shannon entropy for EEWPS distribution is obtained as 
H.M) = - log(a/3A) - ^E[C'(A(1 - £-"•''«>)«) log(C'(A(l - e-ff(5':e))«))] 
+ log[C(A)] -[13- 1)^®[C'(A(1 - iog(i _ 


A 


C(A) 


E[C'(A(1 - log(/i(y; 0))]. 


(4.6) 


Also, the mean residual life function of X is given by 

OO ^ 

r+oo, .Af/ AJ C{X) PnJt °° zg3n)iz]a,nl3,@)dz 

Jt (x - t)f{x)dx _ ^ ^ 


m{t) = E[X - t| A >t] = 


1 - E{t) 

OO 

C(A) E PnE[ZI^z>t)\ 

n=l 


C(A) -C(AG(x)) 


— t 


-1, 


C{X) - C{XG{x)) 
where Z has a EEW distribution with parameters a, n/3 and 0. 


(4.7) 


4.4 Reliability and average lifetime 


In the context of reliability, the stress - strength model describes the life of a component which 
has a random strength X subjected to a random stress Y. The component fails at the instant 
that the stress applied to it exceeds the strength, and the component will function satisfactorily 
whenever X > Y. Hence, R = P{X > y) is a measure of component reliability. It has many 
applications especially in engineering concept. Here, we obtain the form for the reliability R 
when X and Y are independent random variables having the same EEWPS distribution. The 
quantity R can be expressed as 


R = 


roo POO 

/ f{x] a, f3, X,@)E{x‘,a,/3, X,Q)dx = / Xg{x) 
Jo Jo 

OO „QQ 

'^Pn 9{n){x]a,nl3,e 

Jo 


C'{XGix))C{XGix)) 

C'2(A) 


dx 


e^ ciXGjx)) 

U) dx. 


(4.8) 



















5 Estimation 


In this section, we first study the maximum likelihood estimations (MLE’s) of the parameters. 
Then, we propose an Expectation-Maximization (EM) algorithm to estimate the parameters. 


5.1 The MLE’s 


Let xi,..., Xn be observed value from the EEWPS distribution with parameters ^ = (a, /3, A, 
0)^. The log-likelihood function is given by 

n 

ln = ln{^\x) = n[log{a)+log{(3)+ log{X)-log{C{X))] + '^log[h{xi;Q)] 

i=l 

n n n 

-a ^ H{xi-, 0) -h (/3 - 1) ^ log ti + Y^ log(C"(Atf)), 

2 = 1 2 = 1 2=1 


where x = (xi,..., Xn) and ti = 1 — e ^he components of the score function U{^; x) = 


( din 
1 da 


dl^ 

dp 


din 

d\ 


din ' 

< 9 ©- 


are 


din 

da 

din 

d/3 

din 

dX 

din 

dOk 


n 

a 


^iL(xi;0), 

2=1 


'^ + J2^og{ti) + Y 


Atf log{ti)C"{Xtf) 
C'(Atf) 


n nC'jX) ^ t^C'jXt^) 

A C(A) ^ C'(Atf) 


^9/i(xi;0) 1 ^9Lf(xj;0) 

^ dQk h{xi]®) ^ dOk 


+(/3-1)E 

2=1 


dtj 

dQk 

ti 


+ /3^E 




2=1 


c'{xi 


(5.1) 

(5.2) 

(5.3) 


(5.4) 


where 0fc is the A:th element of the vector 0. 

The MLE of 4, say is obtained by solving the nonlinear system U{^;x) = 0. We cannot 
get an explicit form for this nonlinear system of equations and they can be calculated by using 
a numerical method, like the Newton method or the bisection method. Only, for given 0, from 
(j5.1|) we have 

n 

“ " EILi 0) ■ 

Therefore, (|5.4I) becomes 
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dh{xi]®) 1 
dQk 'h{xi;&) 


n 

Eti ®) 
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dH{xi;@) 

dQk 
















(5.5) 


dtj 

de. 


+(/3_l)^i^ + /3A^ 

T'i 




i.i - ..1 ^{Atf) 

Theorem 5.1. The pdf, /(x|©), of EEWPS distribution satisfies on the regularity condistions, 

i.e. 

i. the support of f{x\Q) does not depend on ©, 
a. /(x|©) is twice continuously differentiable with respect to ©, 

Hi. the differentiation and integration are interchangeable in the sense that 

a roo i-oo a q2 roo poo o2 


d&de^- 


Proof. The proof is obvious and for more details, see 


Casella and Bergeii (j200lh Section 10. □ 


The asymptotic confidence intervals of these parameters will be derived based on Fisher 
information matrix. It is well-known that under regularity conditions, the asymptotic distribu¬ 
tion of y/n{^ — is multivariate normal with mean 0 and variance-covariance matrix 
where Jn(^) = lim„_>.oo and In{$) is the observed information matrix as 


In {$) = - 


Uaa UaX 

UaP Up0 Upx 
Uap Uxp Uxx 


T n 


Ujs 

uL 


UaO Up® Ux& I U&& 

whose elements are obtained by derivative the equations (I5.1l) - (|5.4p with respect to parameters. 


5.2 EM-algorithm 

The traditional methods to obtain the MLE’s are numerical methods for solving the equations 
(153])-(lEl, and sensitive to the initial values. Therefore, we develop an EM algorithm for 
obtaining the MLE’s of the parameters of EE WPS class of di s tribu tions. It is a very powerful 


tool in handling the incomplete data problem (jPempster et al 


19771 1. It is an iterative method. 


and there are two steps in each iteration: Expectation step or the E-step and the Maximization 
step or the M-step. The EM algorithm is especially useful if the complete data set is easy to 
analyze. 

Using (|2.3|) . we define a hypothetical complete-data distribution with a joint pdf in the 
form 

a 

g{x, z] $) = ^ za/3h{x; ©)(1 — t)ffl~^, x > 0, 2 G N, 

C(A) 














where t = 1 — The E-step of an EM cycle requires the expectation of 

where 0^^^) is the current estimate (in the rth iteration) of The 

expected value of Z|X = x is 


E{Z\X = x) 


At“C"(At^) 


(5.6) 


The M-step of EM cycle is completed by using the MLE over 0, with the missing z'’s 
replaced by their conditional expectations given above. Therefore, the log-likelihood for the 
complete-data y = (xi,..., Xn, zi ,..., Zn) is 


n n n 

l*{y;^) oc log(A)nlog(a/3)^log/i(xi; 0) - Q:^i?(xi; 0) 

2 = 1 2 = 1 2=1 

n 

+ - 1) log(l - - nlog(C(A)). (5.7) 

2 = 1 


On differentiation of (15.7p with respect to parameters a, /3, A and 0^, we obtain the components 
of the score function as 


2=1 2=1 

It = ^ + 

dll 


^ A C(A) ’ 
dh{xi]®) 1 


E 


-°E ie +Efe<^-‘)- 

i=l i=l 


dQk ^ dQk 'h{xi;G) 50^ ’ i _ e-«//(x,;©)' 


Therefore, we obtain the iterative procedure of the EM-algorithm as 


/3(i+i) 


A(i+i) 


—n 
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dh{xi] 0 


(i+i)x 
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2=1 
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577(xi;0^^ 0 - 1 


2=1 
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'X _ e-aO)//(xi;©^"+^^) 


= 0 , 




















where and 0^ are found numerically. Here, we have 


iy^ = l+ _S i = l,2,...,n, 




where — e 


0 ■)^ 


We can use the results of 


LouisI (Il982l l to obtain the standard errors of the estimators from 


the EM-algorithm. Consider = E{Ic{&',y)\x), where Ic{&',y) = — 


_ fdUiy,&)^ 


a© 


is the 


(/c + 3) X (/c + 3) observed information matrix. If £m(0; x) = Var[U{y\ 0)|aj], then, we obtain 
the observed information as 


J{&\x) = 4 ( 0 ;®) - -^^( 0 ;®)- 


(5.8) 


The standard errors of the MLE’s based on the EM-algorithm are the square root of the 
diagonal elements of the J(0; ®). The computation of th ese matrices are too long and tedious. 


Therefore, we did not present the details. Reader can see 
calculate these values. 
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6 A real example 


In this section, we analyze the real data set given by iMurthv et al.l (|2004l l to demonstrate the 
performance of EEWPS class of distributions in practice. This data s e t cons ists of the failure 


times of 20 mechanical components, and is also studied by I Silva et al 


(120131 1: 


0.067, 0.068, 0.076, 0.081, 0.084, 0.085, 0.085, 0.086, 0.089, 0.098 
0.098, 0.114, 0.114, 0.115, 0.121, 0.125, 0.131, 0.149, 0.160, 0.485 

Since the EEWPS distribution can be used for modeling of failure times, we consider this 
distribution for htting these data. But, this distribution is a large class of distributions. Here, 
we consider five sub-models of EEWPS distribution. Some of them are suggested in literature. 

i. The exponentiated Weibull geometric (EWG) distribution, i.e. the EEWPS distribution 
with H{x, 0) = x'^ and C'(A) = A(1 — A)“^. 

ii. The complementary Weibull geometric (CWG) distribution, i.e. the EEWPS distribu- 
tion with H(x,&) = C(A) = A(1 — A)“^ and /? = 1. This distribution is considered by 
Cordeiro and Silva (2014). 

hi. The generalized exponential geometric (GEG) distribution, i.e. the EEWPS distribution 


with H (x, Q) = X and C'(A) = A(l—A) This distribution is considered by 

(1201J). 


Mahmoudi and Jafari 

























iv. The exponentiated Chen logarithmic (ECL) distribution, i.e. the EEWPS distribution with 
H{x, 0) = exp(x^) and C'(A) = — log(l — A). 

iv. The complementary Chen logarithmic (CCL) distribution, i.e. the EEWPS distribution 
with H(x,Q) = expfx'^), C'(A) = — log(l — A) and /3 = 1. This distribution is considered by 


Cordeiro and Silval (j2014l i. 


The MLE’s of the parameters for the distributions are obtained by the EM algorithm 
given in Section [5l Also, the standard errors of MLE’s are computed and given in paracenteses. 
To test the goodness-of-fit of the distributions, we calculated the maximized log-likelihood 
(log(L)), the Kolmogorov-Smirnov (K-S) statistic with its respective p-value, the AIC (Akaike 
Information Criterion), AICC (AIC with correction), BIC (Bayesian Information Criterion), 
CM (Cramer-von Mises statist ic) and AD (Anderson-Darling sta tistic) for the five submodels of 


distribution. The R software (jR Development Core Teaml . 


201 4 ) is used for the computations. 


The results are given in Tabled and from K-S, it can be concluded that all five models 
are appropriate for this data set. But, the EWG and ECL distributions are better than other 
distributions. In fact, we have a better fit when there is the parameter j3 (exponentiated 
parameter) in model. The plots of the densities (together with the data histogram) and cdf’s 
given in Eigure[T] confirm this conclusion. 


Table 2: Parameter estimates (standard errors), K-S statistic, p-value, AIC, AICC, BIC, CM 
and AD for the data set._ 


Distribution 

EWG 

CWG 

GEG 

ECL 

CCL 

a 

28.665 

25.972 

27.752 

17.111 

22.019 

(s.e.) 

(4.617) 

(11.093) 

(6.841) 

(1.572) 

(10.177) 

7 

0.199 

1.642 

— 

0.136 

1.586 

(s.e.) 

(0.052) 

(0.407) 

— 

(0.026) 

(0.231) 

A 

0.136 

0.012 

0.001 

0. 146 

0.261 

(s.e.) 

(0.918) 

(1.122) 

(0.658) 

(0.032) 

(0.332) 

P 

5.5 

— 

13.825 

7.4 

— 

(s.e.) 

(1.6 e®) 

— 

(8.471) 

2.3 

— 

log (L) 

37.978 

26.422 

32.976 

37.794 

25.759 

K-S 

0.124 

0.264 

0.160 

0.121 

0.262 

p-value 

0.917 

0.122 

0.683 

0.931 

0.127 

AIC 

-67.957 

-46.845 

-59.952 

-67.588 

-45.518 

AICC 

-65.29 

-45.345 

-58.452 

-64.922 

-44.018 

BIC 

-63.974 

-43.858 

-56.965 

-63.606 

-42.531 

CM 

0.048 

0.436 

0.153 

0.051 

0.463 

AD 

0.402 

2.537 

1.136 

0.423 

2.663 



















Histogram 


Empirical Distribution 




X X 

Figure 1: The histogram of the data set with the estimated pdf’s (left), the empirical cdf of 


the data set, and estimated cdf’s (right) for fitted of five submodels. 
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